function [alfa_1 alfa_2 beta_1 beta_2] = func_alfa_beta_from_nz(nz, CONSTS)

    [nxy1 nxy2] = func_nxyk_from_nz(nz, CONSTS); 
    [n1 n2] = func_nk_from_nz(nz, CONSTS);

    eps = CONSTS.eps;
    g = CONSTS.g;
    
    alfa_1 = (nz.^2 + nxy1.^2 - eps)./g;
    
    alfa_2 = (nz.^2 + nxy2.^2 - eps)./g;
    
    beta_1 = (nz./n1);
    
    beta_2 = (nz./n2);
    
end
